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I.  INTRODUCTION 


During  the  past  decade,  there  have  been  many  advancements  in  the  de¬ 
velopment  of  highly  effective  tactical  weapons  for  use  at  intermediate 
range  and  with  high  average  flight  speed.  One  of  the  simplest  is  the 
solid  fuel  ramjet  (SFRJ). 

The  solid  fuel  ramjet  most  often  consists  of  a  solid  fuel  grain  which 
provides  the  walls  for  the  combustion  chamber  [Ref.  1].  A  sudden  expan¬ 
sion  at  the  air  inlet  end  of  the  combustor  can  be  used  to  provide  flame 
stabilization  in  the  SFRJ  (Fig.  1).  Combustion  can  be  sustained  through¬ 
out  the  fuel  grain  provided  that  the  step  height  (h)  is  sufficiently 
large.  However,  the  larger  is  the  step  height,  the  larger  the  losses  in 
inlet  stagnation  pressure. 

Due  to  the  sudden  expansion,  two  distinct  flow  fields  are  generated 
within  the  fuel  grain.  In  the  first  one,  the  recirculation  zone,  the  flow 
in  highly  turbulent  and  normally  fuel  rich.  This  hot  gas  region  provides 
the  energy  necessary  for  sustaining  the  combustion  process  which  must 
occur  further  downstream.  Downstream  of  the  flow  reattachment  a  turbu¬ 
lent  boundary  layer  develops  and  the  combustion  is  similar  to  that  in  a 
hybrid  rocket,  namely  diffusion  controlled.  A  diffusion  flame  emanates 
from  the  recirculation  zone  and  remains  within  the  developing  turbulent 
boundary  layer  between  the  fuel  rich  zone  near  the  wall  and  the  oxygen 
rich  central  core.  Due  to  that  diffusion  flame,  heat  is  transferred  by 
convection  and  radiation  to  the  solid  surface  which  causes  vaporization 
of  the  fuel.  The  radiative  and  convective  heat  transfer  mechanisms  are 
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not  independent.  The  presence  of  radiation  tends  to  increase  the  rate  of 
vaporization  and  this  in  turn  tends  to  reduce  the  heat  transfer  to  the 
wall  by  convection. 

Some  of  the  fuel  leaving  the  grain  surface  fails  to  reach  the  flame 
by  the  end  of  the  fuel  grain,  remaining  instead  in  the  boundary  layer 
zone  below  the  flame.  This  stratified  boundary  layer  exits  from  the  aft 
end  of  the  fuel  grain  and  may  result  in  decreased  combustion  efficiency  if 
the  fuel  is  not  consumed  before  entering  the  exhaust  nozzle.  In  order  to 
increase  the  combustion  efficiency,  an  aft  mixing  chamber  downstream  of 
the  fuel  grain  (Fig.  1)  is  often  employed.  This  allows  additional  chemi¬ 
cal  reaction  of  the  unburned  fuel  to  take  place.  To  promote  mixing  in 
that  region  and  to  increase  the  allowable  fuel  loading,  a  portion  of  the 
inlet  air  is  sometimes  bypassed  to  the  region  aft  of  the  fuel  grain. 

There  have  been  continuing  research  efforts  concerned  with  the  com¬ 
bustion  behavior  of  the  SFRJ  at  the  Naval  Postgraduate  School ,  the  Naval 
Weapons  Center,  and  the  United  Technologies  Corporation/CSD.  At  NPS, 
both  mathematical  modeling  and  experimental  efforts  [Ref.  2,  3,  4,  5,  6] 
have  been  conducted  to  determine  the  effects  of  design  and  operational 
variables  on  the  obtainable  performance. 

The  computer  model  simulation  of  the  SFRJ  combustion  process  has 
evolved  from  an  original  stream- function  vorticity  formulation  [Ref.  1] 
to  a  primitive-variable  (pressure,  velocity)  model  which  includes  the 
aft  mixing  chamber  [Ref.  7],  These  models  have  not  included  radiative 
heat  transfer  to  the  fuel  surface.  This  has  limited  the  utility  of  the 
model  since  many  all -hydrocarbon  fuels  produce  significant  amounts  of 
radiative  heat  transfer  through  the  generation  of  carbon  particulates  in 
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the  flame  zone.  The  purpose  of  the  present  effort  was  to  incorporate 
into  the  primitive  variable  computer  model  radiative  heat  transfer  from 
carbon  particles  within  the  flame  zone  to  the  fuel  surface. 
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II.  DESCRIPTION  OF  COMPUTER  MODEL 


A.  INTRODUCTION 

At  the  Naval  Postgraduate  School  there  have  been  two  basic  computer 
codes  used  to  model  the  flow  within  solid  fuel  ramjets.  The  first  one 
was  based  on  vorticity  (w)  and  stream  function  (<fi)  [Ref.  1].  It  pre¬ 
dicted  reasonably  accurate  velocity  and  temperature  distributions  but 
did  not  predict  pressure  distributions  accurately  and  was  not  easily 
extended  to  complex  geometries  or  to  three  dimensional  flows.  The  second 
computer  model  is  based  on  the  primitive  variables,  velocity  (u,v)  and 
pressure  (p)  [Ref.  7].  It  has  also  been  used  to  predict  the  flowfield 
within  the  fuel  grain.  However  the  u-v-p  model  also  permitted  the  pre¬ 
diction  of  the  flowfield  within  the  aft  mixing  region.  The  latter  model 
can  be  used  to  predict  the  effects  of  fuel  properties  on  the  SFRJ  per¬ 
formance  and  to  evaluate  the  effects  of  different  geometries  and  operating 
conditions.  However,  only  convective  heat  transfer  to  the  fuel  surface 
has  been  considered  and,  therefore,  the  predicted  fuel  regression  profiles 
have  not  been  correct  for  fuel  systems  (such  as  HTPB)  that  generate  large 
amounts  of  radiative  heat  transfer  to  the  fuel  surface. 

In  the  present  work,  the  primitive  variable  computer  program  was  modi¬ 
fied  to  include  the  effect  of  radiation  to  the  fuel  surface  from  carbon 
particles  within  the  flame  zone.  Gas  phase  radiation  was  neglected. 


B.  OVERVIEW  OF  PRIMITIVE  VARIABLE  COMPUTER  PROGRAM 
1 .  Assumptions 

The  CHAMPION  2/E/FIX  computer  program  developed  by  Pun  and 
Spalding  [Ref.  8]  was  used  as  the  basis  for  the  primitive  variable 
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model.  The  flow  was  assumed  to  be  steady,  subsonic,  recirculating, 
axisymmetric  and  to  have  constant  specific  heat.  To  calculate  the 
effective  viscosity,  a  modified  Jones-Launder  [Refs.  7  -  10]  two- 
parameter  turbulence  model  was  used.  The  combustion  was  considered  to 
be  mixing  limited  (infinitely  fast  kinetics)  with  a  simple  one-steD 
chemical  reaction  of  the  form: 

1  kg  fuel  +  i  kg  Oxidizer - ►  (1  +  i)  kg  Products 

Four  species  were  considered:  oxygen,  nitrogen,  fuel  and  products.  Be¬ 
cause  of  the  infinitely  fast  kinetics  assumed,  no  oxygen  could  exist  at 
the  fuel  surface  and  that  surface  was  considered  to  be  isothermal.  The 
turbulent  Lewis  number  was  also  taken  to  be  unity. 

2.  Governing  Equations  and  Boundary  Conditions 

With  the  above  assumptions,  the  governing  equations  for  axisym- 
metrical  flows  can  be  cast  in  the  following  format: 


h  +  F57  <»"•>  -  h  fr.  W  -  F5F  !?!  *  s,  °  0  <’> 


Convection  Terms  Diffusion  Terms  Source  Terms 

m 

where  *  stands  for  the  variables  u.v.k.e.h.x®  (i”fu  -  and  mN2 
($  =  1  for  the  continuity  equation).  S#  is  the  "source  term"  and  r#  is 
the  effective  exchange  coefficient  for  turbulent  flow  (Table  I). 


ox 


Since  the  quantities  h  (stagnation  enthalpy),  mN  and  x  =  j 

have  identical  governing  differential  equations  and,  in  appropriate  dimen¬ 
sionless  form.  Identical  boundary  conditions,  only  one  of  the  equations 
had  to  be  solved. 


verning  Equation  Parameters 


The  dimensionless  enthalpy  was  selected  to  be: 


H  = 


Nn~h 

hin  *  *fg 


(2) 


where  h^  is  the  enthalpy  deep  in  the  fuel  grain  and  h^  the  enthalpy  of 
the  air  at  the  inlet  dump  plane.  The  corresponding  dimensionless  forms  for 
the  other  properties  were: 


m, 


%  -  V 

*in  c 

N2  mN~  -  V 
*in  ^fg 


(3) 


and: 


(4) 


in 


The  solution  of  the  governing  differential  equation  for  h  permitted  calcu¬ 
lation  of  H  which  was  numerically  equal  to  mM  and  x-  The  values  of  mM 

n2  n2 
and  X  were  then  obtained  using  equations  (3)  and  (4). 


From 


2  2 

h  =  h  +  - — tr^—  +  k 


(5) 


h  was  determined.  This,  in  turn, permitted  the  calculation  of  temperature 
(T)from  the  relation: 


mnxAH 

h=  -nr-*'.  (T-Tref> 


(6) 


Equation  (6)  is  a  simplified  method  for  accounting  for  the  heats  of  for- 

mation  and  results  when  all  of  the  heat  of  combustion  (J/._  )  is 

Kgfuel 

"attached"  to  the  oxidant  [Ref.  11]. 
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The  density  was  calculated  from  the  perfect  gas  law: 


(7) 


in  which  R  was  calculated  using  the  average  molecular  weight  (R  * 


I) 


The  mass  fractions  of  fuel  (m^),  oxygen  (mQX)  and  products  (mpr) 


were  calculated  from  the  relations: 

X  >  0;  mfu 

X  <  0;  mfu 

and:  m _  =  1  -  m 

pr 


*  X  .  mox  ■  0 

*  0.  %  =  -XT 


ox  ‘  mfu  ‘  mN, 


(8) 

(9) 


In  order  to  calculate  the  wall  shear  stress  (x  ),  a  two  part  boundary 

W 

layer  was  used.  The  division  between  the  laminar  sublayer  and  the  turbu¬ 
lent  layer  was  assumed  to  be  at  Yp+  *  11.5  [Ref.  8]. 

Then, if  Y  +  >  11.5,  the  wall  shear  stress  (xj  was  calculated  in  terms 
p  -  w 

of  the  near-wall  kinetic  energy  of  turbulence  using  the  relation: 


T 


W 


(10) 


whereas,  if  yp+  < 
the  formula: 


11.5, Tw  was  found  assuming  a  linear  velocity  profile  from 


Tw  "  u 


lam  5 


(11) 


x  provided  the  boundary  conditions  for  the  u  and  v  equations. 

W  -  m 

Since  h,  mN  and  X  =  mf(J  -  -j2-  are  considered  to  be  "conserved  proper¬ 
ties"  [Ref.  11],  the  boundary  condition  for  all  of  them  can  be  written  in  the 
form: 


8 


(12) 


m 


bw 


(F*  17>bw 
cbw  ‘  cfg 


where  $  stands  for  h,  mN  or  X 
c  r<2 


The  mass  transfer  conductance  (g)  is  generally  defined  to  be: 


3$. 


_  (  $  3r  )bw 


$  -  $ 
c  c 


(13) 


bw 


g  can  also  be  equated  to  -£■  [Ref.  7],  Equations  (12)  and  (13)  result  in: 

Up 


$  $ 

cp~  cbw 
bw  ^  $  $ 

cbw  ‘  cfg 


(14) 


where  the  near-wall  grid  value  (*  )  has  been  used  in  place  of  the  Tree 

P 

stream"  value  (*  ). 

'**oo 

The  mass  transfer  parameter  is  given  by: 


BP  - 


$  $  ~  " 

CP  "  Sw  =  hp  '  hbw 

cbw  -  cfg  ^bw  '  **fg 


(15) 


Therefore  (14)  and  (15)  result  in: 


"bw  1  ’BP 


(16) 


Solution  of  the  stagnation  enthalpy  allowed  calculation  of  the  blowing 
parameter  using  (15). 
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The  wall  shear  stress  calculated  for  no  wall  mass  addition  ((10)  or  (11)) 
was  then  modified  for  blowing  using  the  Couette  flow  approximation  [Ref.  11]: 


=  In  (1  +  BP) 
bw  ~  w  BP 


07) 


g  can  then  be  calculated  from  f^/Up  and  mbw  calculated  using  (16). 

The  convective  flux  (qcon)  on  all  boundaries  was  calculated  using  a  modi¬ 
fied  Reynolds  analogy: 


Rcon 


■  %  -  *w>  if 


08) 


where  the  values  of  h  and  u  at  the  near  wall  node  point  (p)  have  been  used 
in  place  of  the  values  at  the  edge  of  the  boundary  layer. 

The  fuel  regression  rate  (RR)  and  the  blowing  velocity  (Vbw)  were 
calculated  from: 


°fg 


09) 


bw 

bw 


(20) 


3.  Solution  Procedure 

Using  equation  (1)  in  finite  difference  form,  the  variables:  u, 
v,  k,  e  and  h  were  found.  The  line  by  line  iterative  procedure  employed 
upwind  differencing  and  under  relaxation  to  promote  convergence  [Ref.  8]. 
Imposing  mass  and  momentum  conservation  on  each  radial  grid  line  and  in 
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each  nodal  control  volume,  the  pressure  (relative  to  a  selectable  position 
and  magnitude  within  the  grid)  was  obtained.  The  effective  viscosity, 
density  and  temperature  were  also  obtained  as  described  above.  More  de¬ 
tails  of  the  solution  procedure  can  be  found  in  Reference  8. 

C.  MODIFICATION  TO  INCLUDE  RADIATION  TO  THE  FUEL  SURFACE 
1 .  Introduction 

The  regression  rate  of  the  fuel  is  proportional  to  the  total  heat 
transfer  to  the  surface.  This  heat  transfer  consists  of  two  interdependent 
parts,  convection  and  radiation.  The  presence  of  radiation  tends  to  in¬ 
crease  the  rate  of  fuel  decomposition/vaporization  and  this  in  turn  tends 
to  reduce  the  heat  transfer  to  the  wall  by  convection. 

The  convective  heat  transfer  depends  on  the  flow  conditions  within 
the  port  of  the  ftel  grain;  primarily  upon  the  air  temperature  and  the  air 
mass  flow  per  unit  irea. 

The  radiant  energy  to  the  fuel  surface  can  be  approximated  by: 


qrad 


(21) 


where:  e,, 
w 

a 


the  emissivity  and  absorptivity  of  the  wall 

the  Stefan-Boltzmann  constant  (=  5.669  x  10"8  -»*--*■) 

m  k 


Cg  =  the  emissivity  of  the  flame  zone 


Tw  =  the  wall  temperature 
Tf  *  the  flame  temperature 

This  simplified  approach  was  taken  in  keeping  with  the  approaches  taken  for 
the  combustion  process  and  the  wall  mass  addition. 
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Radiative  heat  transfer  is  a  function  of  the  pressure  (p)  in  the 
combustor,  the  thickness  (6)  of  the  combustion  (flame)  zone,  the  amount  of 
soot,  and  is  strongly  dependent  on  the  flame  temperature  • 

2.  Assumptions 

The  emissivity  of  the  flame  zone  can  be  written  [Ref.  12]: 

eg  =  1  -  e“aniS  (22) 

where  a  is  an  empirical  constant,  5  is  the  flame  zone  thickness  and  n  is 
the  number  density  of  particulates  in  the  flame  zone(Fig.  2).  In  this  applica¬ 
tion, only  the  cbminant  effect  of  radiation  from  the  carbon  particles  within  the 
flame  zone  was  considered.  It  was  assumed  that  the  gases  between  the 
flame  and  the  wall  were  transparent  and  that  gas  phase  radiation  was  negli¬ 
gible. 

The  number  density  n  of  particulates  in  the  flame  zone  was  assumed 
to  be  proportional  to  the  number  density  within  the  boundary  layer  and 
this  in  turn  proportional  to  the  mass  density  of  unburned  fuel  in  the 
boundary  layer.  The  latter  is  an  assumption  that  a  fixed  percentage  of 
the  unburned  fuel  is  solids.  Thus: 


nflame  % 


nboundary 

layer 


Mass  of  unburned  fuel  _  fu 


W 


VoT 


(23) 


But: 


wfu  3  wmix'mfu 


(24) 


where  wmix  Is  the  mass  of  the  mixture  and  mfu 
fuel.  Also: 


is  the  mass  fraction  of 


wmix 


p(Vol) 


(25) 
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where  p  is  the  gas  density. 

Substituting  (24)  and  (25)  into  (23)  yields: 

n  %  omfu  (26) 

Equation  (26)  indicates  for  the  simplifications  employed  in  this  applica¬ 
tion  that  the  number  density  of  particulates  (n)  depends  directly  on  the 

D 

pressure  (p  =  ^). 

The  flame  zone  thickness  6  has  been  found  to  be  proportional  to 
the  square  root  of  the  port  diameter,  but  for  a  given  geometry  is  almost 
constant  [Ref.  12].  Therefore,  the  effect  of  6  can  be  incorporated  into 
the  empirical  constant  a.  When  these  approximations  are  incorporated  into 
(22)  the  following  expression  is  obtained: 

.  pmfu 

eg  ^  1  -  e 

wfu/Vol 

or  1  .  e  X  (27) 

where  X  an  empirical  constant. 

To  estimate  the  mass  of  unburned  fuel  in  the  boundary  layer  at  a 
particular  axial  location,  a  summation  was  made  over  control  volumes  of 
the  finite  difference  grid  from  the  fuel  surface  to  the  axis  of  symmetry. 
Thus: 

wfu  =  cellsWfucell  =  cellsWmi'xmfu  =  cellsPALmfu  (28) 

Equation  (23)  was  divided  by  the  summation  of  all  control  volumes  that 
had  finite  amounts  of  fuel  present  and  then  substituted  into  (27)  to  esti¬ 
mate  eg. 
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The  mass  transfer  number  BP  must  also  be  modified  to  account  for 


thermal  radiation.  The  modified  form  of  the  definition  of  BP  for  systems 
involving  radiation  heat  transfer  [Ref.  11  ]  is: 


BP 


hbw  "  ^fg 


-  (qrad) 
v  m" 


(29) 


In  this  investigation  Twa^i  was  taken  to  be  Tref  and,  therefore,  (from 
(5)  and  (6)):  hfaw  =  0. 

Thus: 


BP 


- 


(30) 


.  >  >  .  ^  ' 
where  qrad  is  the  radiative  heat  flux  (21)  and  m  the  wall  mass  flux. 

3.  Solution  Procedure 

The  same  solution  procedure  discussed  above  for  the  simple  con¬ 
vection  model  was  followed  in  determining  the  five  variables  (u,  v,  k,  e, 
h).  The  radiative  heat  flux  was  calculated  from  (21)  for  fixed  values  of 
e  ,  T,  and  T  .  A  corrected  blowing  parameter  was  then  found  from  (30), 

W  T  W 

resulting  in  a  modified  mass  transfer  rate  (eq  (16)).  The  modified  re¬ 
gression  rate  (RR)  and  blowing  velocity  (vfaw)  were  obtained  from  equations 
(19)  and  (20). 

In  order  to  provide  numerical  stability  the  effects  of  radiation 
could  be  omitted  for  a  specified  number  of  initial  sweeps  through  the  grid 
field.  This  allowed  the  convective  heat  transfer  solution  to  converge 
before  introducing  the  radiative  effects. 


III.  EXPERIMENTAL  MEASUREMENT  OF  RADIAL  TEMPERATURE  PROFILES 


Two  experiments  were  performed  in  order  to  obtain  two  different  radial 
temperature  profiles  near  the  aft  end  of  the  fuel  grain  for  two  different 
air  flow  rates.  The  intent  was  to  obtain  experimental  profiles  that  could 
be  compared  with  computer  generated  profiles  that  are  in  part  dependent 
upon  air  flow  rate  as  a  data  input.  The  comparison  could  then  be  used  to 
assist  in  the  validation  of  the  computer  model. 

Two  Plexiglas  fuel  grains  were  identically  prepared  with  eight  chromel- 
alumel  thermocouples  imbedded  0.0254  m  from  the  aft  end  of  the  grain.  Four 
of  the  eight  thermocouples  used  36  gage  wire  while  the  remaining  four  were 
30  gage  (Figure  3). 

The  four  thermocouples  of  each  size  were  imbedded  to  distances  of 
0.51,  2.54,  4.57,  and  6.35  mm  from  the  inside  grain  wall.  The  thermo¬ 
couples  were  also  spaced  45  degrees  apart  from  each  other. 

One  test  was  conducted  with  50%  byDass  air  and  one  without  bypass  air. 
The  total  nominal  air  flow  rate  was  0.0907  kg/s.  The  voltage  signals 
obtained  from  the  thermocouples  were  input  directly  into  a  Visicorder. 
Temperatures  were  obtained  from  the  data  by  comparison  of  the  recorded 
signals  with  previously  obtained  calibration  curves  (Figure  4).  The  cali¬ 
bration  curves  (for  both  gages  of  thermocouples)  were  necessary  since  the 
experiments  were  conducted  without  the  use  of  an  ice  bath  or  other  type 
of  reference  junction.  The  calibration  consisted  of  using  an  oxy-acetyl ene 
torch  as  a  heat  source  for  the  thermocouples  that  were  simultaneously  being 
referenced  to  an  electronic  ice  point.  By  measuring  the  temperature  in¬ 
duced  voltages  on  a  digital  voltmeter  the  deflections  on  the  Visicorder 
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traces  could  be  related  to  the  temperature  induced  voltages.  The  data 
extraction  method  consisted  of  measuring  the  deflection  distance  from 
the  Visicorder  trace,  entering  the  calibration  curve  with  the  distance 
and  obtaining  a  corresponding  voltage  which  was  then  used  with  tables 
that  relate  the  induced  thermocouple  voltage  to  the  actual  temperature 
being  sensed  by  the  thermocouple  (uncorrected  for  radiation  heating). 
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IV.  DISCUSSION  OF  RESULTS 


A.  INTRODUCTION 

In  the  modifications  made  to  the  primitive  variable  computer  pro¬ 
gram  only  the  dominant  effect  of  radiation  from  the  carbon  particles 
within  the  flame  zone  was  considered.  The  radiative  heat  transfer  to 
the  fuel  surface  was  then  approximated  by  equation  (21)  where  the  emis- 
sivity  of  the  flame  zone  (eg)  was  taken  from  the  relation  (27).  The 
effects  of  the  amount  of  unburned  carbon  and  the  combustor  pressure  were 
incorporated  into  this  expression.  For  the  Plexiglas-02  system  (which 
has  only  small  radiation  effects)  a  typical  value  for  cg  is  0.02  and  the 
radiative  heat  flux  to  the  fuel  grain  is  approximately  15%  of  the  con¬ 
vective  flux  [Ref.  13].  Using  these  values  the  empirical  constant  (A) 
in  equation  (27)  was  found  to  be  approximately  3.  For  fuel  systems  that 
generate  larger  amounts  of  radiative  heat  transfer  (for  example  HTPB) 
the  empirical  constant  A  should  be  lower  than  3. 

It  was  also  assumed  that  the  radiation  (drad)  is  absorb^  and  emitted 
at  least  an  infinitesimal  distance  inside  the  solid  surface  and  that  the 
gases  between  the  flame  and  the  wall  were  transparent.  The  region  between 
the  flame  and  the  wall  was  taken  to  be  totally  unaffected  by  the  radia¬ 
tion  [Ref.  11]. 

Once  numerically  stable  computer  model  solutions  were  obtained,  the 
effects  of  pressure,  inlet  velocity  and  increased  qrad  (by  reducing  the 
empirical  constant  A  for  systems  that  produce  large  amounts  of  radia¬ 
tion)  were  examined. 
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It  has  been  found  experimentally  that  increased  combustor  pressure 
increases  the  average  regression  rate.  However,  for  PM  the  effect  of 
pressure  alone  is  not  large  [Ref.  2].  Much  of  the  reported  effects  of 
Pc  in  Ref.  2  appear  to  have  resulted  from  variations  in  G  or  inlet 
geometry.  Small  variation  in  r  with  Pc  at  fixed  G  cannot  be  explained 
in  terms  of  the  convection  heat  transfer  because  it  normally  depends 

m 

only  on  mass  flux  (G  =  J-): 


''con 


g0.8  ,  (.m 


0.8 


0.8 


(32) 


For  fixed  values  of  m  and  Ap  (i.e.,  G),  qCQn  is  independent  of  pressure 
since  the  velocity  varies  inversely  with  pressure.  Pressure  changes 
that  result  from  decreasing  the  nozzle  throat  size  necessarily  increase 
qrad*  Therefore»  dra(j  is  thought  to  be  the  cause  for  the  observed  in¬ 
creased  average  regression  rate  with  combustor  pressure.  In  some  com¬ 
bustion  systems  finite  rate  kinetics  can  also  result  in  significant 
effects  of  pressure. 

Fuels  that  generate  large  amounts  of  carbon  (like  HTPB)  also  result 
in  a  different  regression  rate  pattern  compared  with  the  pattern  obtained 
using  PMM.  The  main  difference  is  that  in  the  HTPB  type  systems  the  re¬ 
gression  rates  increase  with  axial  distance  at  the  aft  end  of  the  grain. 
One  possible  cause  of  such  an  effect  is  that  as  the  reacting  flow  gets 
hotter  and  accelerates,  convective  heat  transfer  is  increased.  However, 
the  computer  calculations  Indicated  that  this  was  not  a  dominant  effect. 
Therefore,  q„a(j  was  again  considered  to  be  the  reason  for  this  behavior. 


) 
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The  experimentally  obtained  regression  rate  pattern  for  PfW  [Ref.  7] 
is  higher  than  that  predicted  using  the  primitive  variable  computer  model 


with  only  convective  heat  transfer,  especially  at  the  back  of  the  fuel 
grain.  The  flow  reattachment  point  has  also  been  found  to  occur  slightly 
further  downstream  than  predicted. 

B.  EFFECT  OF  RADIATION  ON  REGRESSION  RATE 

Figure  5  shows  the  regression  profiles  obtained  for  the  two  cases: 
with  and  without  radiation  heat  transfer  to  the  fuel  surface.  The  flow 
reattachment  point  was  not  affected  by  incorporation  of  radiation  into 
the  model.  Radiation  increased  the  regression  rates  downstream  of  re¬ 
attachment,  resulting  in  a  regression  profile  in  better  agreement  with 
experiment. 

C.  EFFECT  OF  PRESSURE 

Figure  6  shows  the  effects  of  combustion  pressure  on  the  fuel  re¬ 
gression  rate  with  radiation.  In  the  first  case  (nominal  Pc)  the 

5  N 

pressure  was  taken  to  be  4.826  x  10  —*  whereas  in  the  other  case  the 

5  N  m 

chamber  pressure  was  9.652  x  10  It  can  be  seen  that  increased 

nr 

chamber  pressure  increased  the  regression  rate  in  general  agreement  with 
experimental  observations.  The  reattachment  point  also  moved  slightly 
to  the  right  in  this  case.  Specific  experimental  data  were  not  availa¬ 
ble  for  comparison. 

As  discussed  above,  the  regression  rate  of  Plexiglas  has  been  re¬ 
ported  to  vary  significantly  with  chamber  pressure  [Ref.  2].  Closer 
examination  of  that  data  reveals  that  for  fixed  geometry  and  mass  flux 


the  pressure  effect  is  quite  weak.  For  the  two  cases  considered,  the 
model  predicted  regression  rate  to  vary  as  Pc^'0^.  Therefore,  the  pre¬ 
dicted  dependence  of  regression  rate  on  the  pressure  was  very  weak  and 
in  reasonable  agreement  with  experiment. 

D.  EFFECT  OF  INLET  AIR  MASS  FLUX  (VELOCITY) 

The  effect  of  increasing  air  mass  flux  on  the  regression  rate  is  shown 
in  Figure  7.  The  regression  profile  remained  the  same  and,  as  expected, 
increased  with  increasing  inlet  velocity,  V^n  (G  ^  V-n).  It  has  been 
found  experimentally  that  the  regression  rate  of  Plexiglas  varies  with 
the  air  mass  flux  (G)  as: 

r  =  C  Gy  (33) 

For  the  three  cases  considered  here  (V.n  =  158,  197  and  236  ^~_)  the 
constant  y  was  predicted  to  be  between  0.23  and  0.28  while  Mady  et  al 
[Ref.  4]  experimentally  found  it  to  be  approximately  0.38.  Therefore, 
the  model  appears  to  correctly  predict  the  nature  of  the  change  in  con¬ 
vective  heat  flux  to  the  fuel  surface  with  air  flow  rate. 

E.  EFFECT  OF  INCREASING  RADIATION 

In  fuel  systems  that  produce  larger  amounts  of  carbon  particulates 
than  the  Plexiglas,  the  empirical  constant  x  must  be  reduced  to  less  than 
3.  For  the  three  test  cases  of  this  study  (Figure  8)  the  values  of  x 
considered  were:  3.00,  2.18  and  1.63.  Lower  values  of  x  resulted  in 
higher  regression  rates  near  the  aft  end  of  the  fuel  grain.  This  results 
from  the  Increasing  amounts  of  unburned  fuel  below  the  flame  with  increasing 
axial  distance. 
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F.  DISTRIBUTION  OF  RADIATION  ALONG  THE  FUEL  GRAIN 

Figure  9  shows  the  predicted  fraction  of  the  radiative  and  convective 
heat  fluxes  to  the  fuel  surface  along  the  grain  for  PW.  From  the  rela¬ 
tions  (30),  (17)  and  (18)  it  can  be  seen  that  the  energy  of  radiation 
tends  to  increase  the  mass  transfer  number  BP,  which  decreases  the  wall 
shear  stress  tl .  This  in  turn  decreases  the  rate  of  convective  heat 
transfer  to  the  wall.  Thus,  increasing  radiation  produces  a  decrease  in 
the  convection  heat  trnasfer  to  the  wall  as  is  depicted  in  Figure  9.  The 
amount  of  radiation  in  the  model  depends  strongly  upon  the  empirical  con¬ 
stant  x.  This  value  is  selected  to  match  experimental  data  for  one  fuel 
and  one  set  of  test  conditions.  More  data  are  required  to  determine 
whether  or  not  one  value  of  X  can  be  used  for  varying  test  conditions 
with  one  fuel. 

G.  COMPARISON  OF  PREDICTED  AND  MEASURED  TEMPERATURE  PROFILES 

1 .  Nonbypass  Conditions 

Figure  10  is  a  graph  of  temperature  vs.  time  for  seven  of  the 
eight  thermocouples  used  in  the  experiment  without  bypass  air  (one  small 
thermocouple  showed  no  response).  The  similarities  in  the  response  times 
of  each  of  the  thermocouples  around  the  circumference  indicated  that  a 
rather  uniform,  ax i symmetric  flame  boundary  existed  at  this  particular 
location  within  the  fuel  grain. 

The  actual  average  regression  rate  for  this  experiment  was  deter¬ 
mined  to  be  0.160  mm/sec,  based  on  the  measured  weight  loss  and  previous 
experiments  using  the  same  configuration.  Using  this  regression  rate  and 
the  information  from  Figure  10,  a  temperature  vs.  distance  profile  was 
produced  as  shown  in  Figure  11.  This  graph  shows  the  similarities  in 
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each  of  the  profiles  for  the  four  heavy  gage  thermocouples.  The  graph 
represents  the  temperature  profiles  as  they  occurred  in  the  fuel  grain 
if  the  ordinate  of  the  graph  is  located  as  shown.  This  would  seemingly 
contradict,  however,  the  knowledge  that  the  surface  temperature  of  burning 
PMM  is  approximately  600  degrees  Kelvin.  In  order  to  take  this  into 
account  the  profiles  need  to  be  shifted  either  to  the  left  or  to  the 
right  within  a  range  of  0.76  mm  depending  upon  the  thermocouple.  Several 
experimental  conditions  affected  the  profiles.  The  leveling  off  of  the 
temperature  occurred  because  the  thermocouples  bent  over  when  they  lost 
strength  as  they  approached  the  flame.  Each  thermocouple  also  generated 
near-wall  mixing  and  increased  the  local  fuel  regression  rate.  This,  in 
turn,  could  reduce  the  local  temperature  because  of  the  more  fuel  rich 
condition  that  was  generated. 

2.  Bypass  Air 

This  test  was  made  so  that  a  lower  air  flow  velocity  in  the  fuel 
port  could  be  achieved  while  still  maintaining  a  choked  exhaust  nozzle. 
However,  bypass  air  may  affect  the  temperature  profile  within  the  fuel 
port.  To  minimize  this  effect  the  motor  air  inlets  were  operated  in  a 
Choked  condition. 

Figure  12  shows  the  temperature  vs.  time  profiles  for  six  of  the 
eight  thermocouples  in  this  experiment.  Two  of  the  thermocouples  were 
never  exposed  to  the  flame  due  to  the  lower  regression  rate  of  the  bypass 
configuration.  Again  there  were  similarities  between  all  of  the  profiles. 

For  this  experiment  the  average  regression  rate  was  determined  to 
be  0.109  mm/sec  based  on  the  measured  weight  loss  and  previous  experiments 
using  the  same  configuration.  Figure  13  is  a  plot  of  temperature  vs. 
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distance  based  on  this  regression  rate.  This  graph  shows  good  agreement 
between  the  profiles  obtained  with  the  heavy  gage  thermocouples.  In  order 
to  bring  these  profiles  to  a  position  that  would  agree  with  a  600  degree 
Kelvin  surface  temperature,  a  shift  to  the  left  of  from  0.57  to  0.70  mn 
in  the  profile  would  be  required. 

3.  Comparison  with  Computer  Model  Temperature  Profiles 

Using  the  primary  air  mass  flow  rate,  the  chamber  pressure  and 
the  temperature  of  the  inlet  air  as  inputs  to  the  computer  program,  tem¬ 
perature  vs.  distance  profiles  were  predicted  as  shown  in  Figures  11  and 
13  by  the  dashed  lines.  The  solid  lines  represent  experimental  profiles. 

The  dashed  line  in  Figure  13  is  based  on  a  primary  air  mass  flow 
rate  of  0.0445  kg/sec.  For  the  range  of  temperatures  between  600  and 
1 000°K  the  profiles  were  similar.  However,  above  1000°K  there  was  a 
marked  difference  between  the  profiles.  This  difference  was  probably 
due  to  the  fact  that  the  thermocouples  had  by  this  time  bent  over  and/or 
burned  out,  making  them  incapable  of  measuring  the  higher  temperatures 
encountered.  Maximum  expected  temperatures  are  approximately  2500°K. 

The  dashed  line  in  Figure  11  was  based  on  an  air  mass  flow  rate 
of  0.088  kg/ sec.  There  was  even  less  agreement  between  profiles  than 
obtained  with  bypass  air  (lower  air  mass  flow  rate).  Again,  this  could 
be  attributed  to  the  fact  that  the  thermocouples  had  ceased  to  function 
at  a  low  temperature  (approximately  1500°K). 

Comparison  between  the  theoretical  and  experimental  profiles  indi¬ 
cates  that  the  computer  model  predicts  a  steeper  gradient  for  both  the 
high  and  low  air  flow  rates  than  those  that  actually  occurred  in  the  ex¬ 
periment. 


The  computer  model  also  predicts  a  steeper  gradient  for  the 
higher  flow  rate,  but  the  experimental  results  indicate  nearly  identical 
gradients  for  both  flow  rates. 

If  the  experimental  profiles  near  the  fuel  surface  are  even 
reasonably  accurate,  it  indicates  that  the  computer  model  predicts  more 
near  wall  mixing  (higher  turbulence)  than  actually  occurs. 
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V.  CONCLUSIONS 


Addition  of  radiative  heat  transfer  to  the  primitive  variable  com¬ 
puter  model  resulted  in  improved  agreement  with  experimentally  obtained 
fuel  regression  rate  profiles.  The  effects  of  varying  chamber  pressure 
and  air  mass  flux  appear  to  be  qualitatively  correct  but  additional  data 
are  required  to  further  validate  the  model. 

Comparison  with  limited  experimental  data  for  radial  temperature  pro¬ 
files  indicates  that  the  model  predicts  too  high  near  wall  turbulence. 
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Fig.  1.  Schematic  of  Solid  Fuel  Ramjet  Combustion  Regions 


SFRJ  motor  g 


BOUNDARY 
LAYER  EDGE 

ESC 

-*4 


^4 


8^ 


'con 


rrad 


DIFFUSION 

FLAME 
n 


777  /  7  /  /7  /  fuel/  ///nr// 


Fig.  2.  Schematic  of  SFRJ  Boundary  Layer  Combustion  Process 
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Figure  3.  Schematic  of  Fuel  Grain  with  Imbedded  Thermocouples. 
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Figure  5.  Fffect  of  Radiation  on  the  Predicted  Fuel  Regression  Rates 
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Figure  6.  Effect  of  Pressure  on  the  Predicted  Fuel  Regression  Rates 
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Figure  7.  Effect  of  Inlet  Velocity  on  the  Predicted  Fuel  Regression  Rates 
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gure  8.  Effect  of  Increasing  Radiation  on  the  Predicted  Fuel  Regression  Rates 
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Figure  9,  Predicted  Distribution  of  Radiation  Along  the  Fuel  Grain 


Figure  10.  Temperature  vs.  Time  (No  Bypass  Air) 


Figure  12.  Temperature  vs.  Time  (Bypass) 
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